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We investigate the entanglement properties of a one dimensional chain of qubits coupled via nearest neighbor 
spin-spin interactions. The entanglement measure used is the ^-concurrence, which is distinct from other mea- 
sures on spin chains such as bipartite entanglement in that it can quantify "global" entanglement across the spin 
chain. Specifically, it computes the overlap of a quantum state with its time-reversed state. As such this measure 
is well suited to study ground states of spin chain Hamiltonians that are intrinsically time reversal symmetric. 
We study the robustness of n-concurrence of ground states when the interaction is subject to a time reversal 
antisymmetric magnetic field perturbation. The n-concurrence in the ground state of the isotropic XX model is 
computed and it is shown that there is a critical magnetic field strength at which the entanglement experiences 
a jump discontinuity from the maximum value to zero. The n-concurrence for thermal mixed states is derived 
and a threshold temperature is computed below which the system has non zero entanglement. 

PACS numbers: 03.65.Ud, 05.50.+q,75.10.Jm 



> 

\o 

O 

o. 

O' 

Ok 

■*-» 

a 
cr. 



* 



I. INTRODUCTION 

There is considerable interest in understanding the distinc- 
tion between quantum and classical correlations in many -body 
systems. Common examples found in nature are collections of 
spins coupled by pairwise interactions. Spin chain Hamilto- 
nians are described by nearest neighbor interactions between 
spins (usually s = 1/2) particles. A typical example is the one 
dimensional quantum XYZ model: 



HxYz(h) = £-/y<x}cx- +1 +J y a y j G 



/+!' 



-J z (fj(fj +l +h(f r (1) 



which describes pairwise interactions with a homogeneous ex- 
ternal magnetic field h that acts locally on each spin. While 
these models are only an approximation to the real physics, 
they are rich enough to extract essential statistical properties 
of the underlying system. A striking phenomenon in many of 
these systems is the existence of a quantum phase transition 
(QPT) described by the analytic discontinuity in some ther- 
modynamic quantity with the variation of a system interaction 
parameter. For the XY model (J z = 0) the QFT is manifested 
by the divergence in the range of pairwise spin correlations at 
a critical magnetic field strength. 

Given that long range classical correlations are present in 
the ground states of spin chains it is to be expected that these 
correlations can be related to functions quantifying quantum 
entanglement. Several significant results have already been 
obtained for the XY model. It has been shown that a measure 
of entanglement between pairs of spins, the 2-concurrence, 
experiences a QFT at the same critical magnetic field strength 
as the transition point for classical spin correlations [jj 0|. 
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These studies also show that in the isotropic, or XX model 
(J x = J y ,J z = 0), the range of pairwise entanglement is infi- 
nite in the thermodynamic limit even though this case does 
not admit a QFT. The entanglement in a bipartite division of 
a spin chain into two contiguous blocks has been studied in 
Ref . y[ ■ There the von Neumann entropy of the reduced state 
of one block was used as the entanglement measure. By this 
measure, that the entanglement obeys universal scaling laws 
in accordance with conformal invariance. In another work, it 
was proven that there is another characterization of entangle- 
ment, the localizable entanglement £,e, whose range is always 
at least as long as classical correlation lengths |4]. The quan- 
tity t,£ is defined as the maximum pairwise entanglement that 
can be localized on two qubits, on average, by optimizing over 
local operations on the other qubits. 

These investigations focused on bipartite entanglement be- 
tween individual spins or blocks of spins. Recently, it was 
shown that the notion of entanglement can be generalized be- 
yond subsystems by computing the purity of state with respect 
to a chosen subalgebra 1 5 ] . Applied to the X Y model, the rele- 
vant subalgebra is the set of number non-conserving fermionic 
operators that connect different irreducible representations of 
the unitary Lie algebra u(n), where n is the number of modes. 
The purity can then be expressed in terms of fluctuations of 
the total fermion number and it is characterized by a second 
order QFT at the critical magnetic field strength |6]. 

All of the aforementioned results add considerable insight 
into how nonlocal correlations between spins are distributed in 
the ground states of spin chain Hamiltonians. In this paper we 
take a different approach to studying the entanglement of spin 
chains. We do not seek a correspondence in the behavior of 
entanglement and classical correlation functions near critical 
points. Rather, we investigate the stability of "global" entan- 
glement across the system when subject to environmentally 
influenced effects such as finite temperature and perturbation 
by a magnetic field. Any approach to this problem must con- 



tend with the non-uniqueness of a single measure of global 
entanglement in a multipartite system. We pick a measure 
of entanglement, the n-concurrence, that is motivated by an 
underlying time reversal symmetry of the XYZ model. The 
hope is that by studying the behavior of global entanglement 
we can gain some insight into how long range quantum cor- 
relations appear in systems that are already highly correlated 
classically. 

The paper is organized as follows. The properties of the 
entanglement measure are discussed in Sec. [II] In Sec.llIIIthe 
ground state entanglement in the one dimensional XY model 
with periodic boundaries is studied. It is shown that the XX 
model displays a jump discontinuity in n-concurrence when 
the perturbing magnetic field reaches a critical strength. The 
entanglement in a thermal state is computed in Sec. II Vi and the 
threshold temperature is derived which sets an upper bound 
on how mixed the state can be before the n-concurrence van- 
ishes. The computation of entanglement for thermal states 
relies on an important theorem derived in Appendix lAl that 
yields a closed form expression for the concurrence of mixed 
states. In Sec. [V]the preceding analysis is extended to the 
quantum XX model with open boundaries. Issues concerning 
the experimental observation of entanglement in spin chains 
are discussed in Sec. IVII Finally, a summary and conclusions 
are presented in Sec. IVIII 



II. PROPERTIES OF ^-CONCURRENCE 

A basic mathematical tool for studying entanglement is the 
entanglement monotone. It is a mathematical function that 
maps states to real numbers and exhibits two important prop- 
erties. First, it is zero for separable states, i.e. those that can 
be described purely by classical probability distributions; sec- 
ond, it is non-increasing on average under local operations 
and classical communication. There are many monotones to 
choose from that quantify entanglement in multipartite sys- 
tems. The choice of measure may be best dictated by the un- 
derlying symmetries of the system if there are any. 

The monotone we study is concurrence. The 2-concurrence 
was originally derived by Wootters [7] and later generalized 
to any even number of qubits n |8]. On a pure state \\\t) the 
n-concurrence is defined 



Cn(W)) = |(¥l«l¥)| 



(2) 



where U is an anti-unitary time-reversal operator. When act- 
ing on n qubits, we can write 13 = [TR =1 (— JCTy)]x, where x 
is the complex conjugation operator. The n-concurrence and 
its square, the n-tangle, have been shown to be entanglement 
monotones |8, 9] for n even, and for n odd are identically zero. 
The range of the measure is < C„ ( |\|/) ) < 1 . 

The n-concurrence is an attractive measure for two rea- 
sons. First, it is sensitive to global entanglement in the 
sense that it is zero if any qubit is separable from the 
rest of the system. This does not imply that each qubit 
is entangled with every other qubit, however. For in- 
stance, the n-party Greenberger-Horne-Zeilinger (GHZ) state 



|GHZ„) = 1/V2(|0...0) + |1...1>) is maximally concur- 
rent but this can also be the case for sub-global entangle- 
ment, e.g. C8(|GHZz|.) (gi IGHZ4)) = 1. Furthermore, some 
entangled states have vanishing concurrence, e.g. \W^) — 
(l/2)[ |0001) + |0010) + |0100) + |1000) ]. Nevertheless, it 
is sensitive to entanglement described by superpositions of 
states and their spin flips. Second, the n-concurrence mea- 
sures the overlap of a quantum state with its time reversed 
state and it is thus a natural function to choose for eigenstates 
of a Hamiltonian that respects this symmetry. Any Hamilto- 
nian that can be written as a sum of tensors of an even number 
of non identity Pauli operators only is time reversal symmet- 
ric. For example, exclusively pairwise interactions satisfy this 
condition 



UH XYZ (0)U- 1 = H XYZ (0). 



(3) 



This symmetry has important consequences with regard to en- 
tanglement of eigenstates. 

Proposition II.l (| 10]) Let H be a Hamiltonian on some 
number n of quantum-bits. Suppose H has time-reversal sym- 
metry with respect to 15. Let X be a fixed eigenvalue of H. 
Then either (i) X is degenerate (i.e. has multiplicity at least 
two,) or (ii) the normalized eigenstate \X) has C n (\X)) = 1. 
Should n = 2p — 1 , p G N, then case (i) always holds. 

Notice that if any spin does not interact with the others in 
a collection of spins then there will be a degeneracy. There 
are several examples of spin-chain Hamiltonians with non- 
degenerate ground states. Among them is the XYZ Hamilto- 
nian with J x = J y =J z = J>0, also known as the Heisenberg 
interaction. It's ground state was proven by Lieb and Mattis 
IU2I1 to be non-degenerate in any number of dimensions, with 
or without periodic boundary conditions provided the under- 
lying lattice has reflection symmetry about some plane. Addi- 
tionally, in ID, the XY Hamiltonians with J x ,J y ^Q,J Z =Q are 
non-degenerate 1 1 3 ] . In this paper we focus on entanglement 
properties of the latter. 

Extending measures of entanglement on pure states to en- 
sembles of pure states, or mixed states, can be carried out by 
averaging over the entanglement of pure states in the ensem- 
ble. The choice of the state decomposition should not in- 
crease the amount of entanglement, therefore, the entangle- 
ment should be minimized over all valid decompositions of 
the state. This formulation is known as the convex roof of the 
function and for n-concurrence it is written: 



C„(p) =min X>;C„(|y»; p = J>ylvWl 



(4) 



Remarkably, the n-concurrence of a mixed state on n even 
qubits can be expressed in closed form (see AppendixlAl: 



f N-\ 

C(p)=max^0,A, - £ Xj 



(5) 



Here {Xj}j~q = spec(-\/pL5pl5 _1 ), where spec(A) denotes 
the spectrum of the operator A, and N = 2" is the dimension of 
the system. The set of eigenvalues are real, positive numbers 
arranged in non-increasing order: Xq > Xi . . . > Xn- 1 . 



III. GROUND STATE ENTANGLEMENT 



The XY model with a uniform magnetic field h is written 



Hxy (h) = £ J x a]<r] + , + J y a]^ + , + ha\ . (6) 

Here we have assumed cyclic boundary conditions (<3„ +k = 
a"). Additionally, we assume n to be even. The pres- 
ence of the magnetic field breaks the time reversal symme- 
try of the interaction. Indeed, the magnetic field interaction 
Hb = h Y."i= i g) contains a sum of terms involving only odd 
numbers of Pauli operators and therefore it is time-reversal 
anti-symmetric 1 10]: UHgU^ 1 = —Hb, i.e. 



UH XY {h)U- 1 = H XY (-h). 



(7) 



We study how robust the n-concurrence of the ground state is 
to perturbation by the magnetic field. 

The eigenvalues and eigenvectors of Hxy (h) can be solved 
for exactly by performing the Jordan-Wigner transformation 
from Pauli operators to fermionic operators followed by a 
Fourier transformation. The Jordan-Wigner transformation is 
given by 



a j= v J°J 



Clj = O+Vj, 



(8) 
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FIG. 1: Plot of the 4-concurrence of the XY model as a function 
of magnetic field strength and anisotropy y= (J x — J y )/2J where 
J = (Jx +-A0/2- In the isotropic case (J x = J y ), there is a critical mag- 
netic field strength h cr n where the n-concurrence experiences a jump 
discontinuity from a value of one for \h\ < h cr it to zero at \h\ > h cr i,. 
For n = 4, h cri , = 27tan(jt/8). 



where 



H k = 



2e k 


-25, 


25, 


2e k 















for 1 < k < n, 



(13) 



where the introduction of the non-local terms V/ 



J-i. 



Cl 



(1 —2n k ) = (— l) L k=i a k a k enforces the anticommutation 



relations: 



{a\,a k i} = h kk >, {a k ,a k >} = {a^ fl !'} = 



(9) 



Define the Fourier transformed creation and annihilation op- 
erators as: 



A\ 



1 



L a k e 



iin(jk/n- 1/4) 



(10) 



k=\ 



where A[ n+k 



\Al. Following Katsura 1 1 3 ] we can decom- 
pose Hxy(h) into two subspaces corresponding to the ± 



^=i u /t 



eigenvalues of the operator v n +i 

(— l)£*=i A 2t A 2* This amounts to a partitioning of the system 
into even and odd parity halves of the total number of excita- 
tions: 

HxY(h)= l -(l+v n+l )H++ l -(\-v n+x )H-. (11) 

In terms of the new fermionic operators A k ,A k , the Hamil- 
tonian Hxy reduces to block diagonal form with each block 
spanned by particle occupation in positive and negative "mo- 
mentum" number states A' k A k . The even parity sector of Hxy 
is 



n/2 

H = 2_ i H 2k -\ 



(12) 



in the basis {|0),A;[A_J0),a£|0),A_J0)}. The matrix ele- 
ments are 



£*: = (Jx + Jy)cos(nk/n) —h, 5, = (J x —J y )sin(nk/n). 



Similarly, the odd parity sector is, 



n/2-l 
H = H Q + H„+ £ H 2k 

k=l 



where the two boundary terms are 
Ho=(- £ ° °), H, 



-en 

En 



(14) 



(15) 



(16) 



k=\ 



in the bases {|0),Aj|0)} and {|0),A^|0)} respectively. 

The Hamiltonian can be brought to full diagonal form by 
a Bogoliubov transformation on the operators H k to a set of 
fermionic operators p,,P| whose number is conserved. The 
transformation is: 

p, = cos(6,.)A A . + sin(e A .)Al A ., pj =cos(9,)AJ; + sin(e,)A_, 

(17) 
where tan(29,) = h k /t k . In many treatments of the statistical 
properties of the XY model, the boundary terms Hq,H„ are 
dropped with the justification that their contribution is negli- 
gible in the thermodynamic limit 1 1 1 ]. Because we are inter- 
ested in entanglement properties of the system for all even n, 
we keep these terms. 

We wish to compute the n-concurrence of the ground state 
of Hxy(h) which we denote \\\t g (h)). One approach is to 



explicitly calculate C n {\\f g (h))) = \(\\f g (h)\l3\\\f g (h))\. Be- 
cause the Hamiltonian is real, the energy eigenstates are 
real and this amounts to computing the expectation value 
( y Vg\(~i® y ) m \ } ¥g}- After performing the transformation from 
the operators o\- to the fermionic operators, there will be 

0(n 2 ) terms to sum in the expectation value. An alterna- 
tive approach is facilitated using the expression for the n- 
concurrence of mixed states, Eq. [5] The problem is reduced 
to finding spec(y / pT5pTS^), associated with the ground state: 



pr=o = \Vg(h))(y g Qi)\ 



lim - 

r->0 



-$H XY (h) 



(18) 



where Z = Tr[<?~P H ] is the partition function. In describing 
the state as the T — > limit of the thermal sample we have 
implicitly assumed the ground state is non-degenerate. All 
ground states of Hxy(0) with J x ,J y ^ are non degenerate but 
degeneracies do occur at finite magnetic field strengths. In the 
present discussion we assume that at any given degeneracy 
point, the state of the system is in a single (zero entropy) pure 
state. This assumption is dropped in Sec. HVI where it is shown 
that thermal states with degenerate ground states have zero n- 
concurrence. Proceeding with this qualification in mind we 
have 



Op^o^T 1 = 15[ lim £ 



-Wxy(K 



O 



-1 






lim 

r->0 

lim 



,-P%[')g' 
z 

e -p%y(-A) 



(19) 



= \v g (-h))(v g (-h)\, 

where in the fourth line we have used Eq. [7] The matrix 
Pr=0^5pr=o^ is at most rank one and the n-concurrence of 
the ground state of Hxy(h) is therefore, 



C»(|Y,(A)» = |<V«(A)IV,H0>| 



(20) 



The entanglement of the ground state of the XY model 
in n = 4 qubits is plotted in Fig. [2 Notice the symme- 
try of the n-concurrence under h — > — h as anticipated in Eq. 
l20l At zero magnetic field Hxy is ti me reversal symmetric 
and non-degenerate so the 4-concurrence is equal to one for 
anisotropies in the range < y < 1 . In the isotropic case, there 
appears to be a jump discontinuity in the entanglement at a 
critical magnetic field strength h cr \ t . We proceed to show that 
this phenomenon arises for all even n and that h C rtt is given 
by the minimum field strength at which the ground eigenstate 
becomes doubly degenerate. 

Henceforth, we focus on the isotropic XX model (J x = J y = 
J). In this case, the Hamiltonian Hxx (n) is already diagonal in 
the number operators AlAjfc,AljA_£, The eigenvalues of H + 
are given by: 



E+ = 



= < £,*2*-i; 



**€{-e*,e*,0,0} 



(21) 



There are a total of 4"' 2 eigenvalues obtained by the sum over 
elements in the set. However, the projection 1/2(1 + v„+i )H + 
requires that each energy eigenvalue be equal to a sum of 
n/2 terms in the set with the difference between the number 
of ± signs in the sum satisfying #(+) — #(— ) = n/2 mod 2. 
This makes the total number of eigenvalues obtained from H + 
equal to 2"/2. Similarly, 



E = I ±e ±E„- 



B/2-1 >| 

- Y, X2k > x k€{-£k,£k,0,0}>, (22) 
A=l j 

where the projection 1/2(1 — v„+i)//~ requires that each en- 
ergy eigenvalue be equal to a sum of n/2 terms satisfying 

# (+) -#(-) = («/2 + 1) mod 2. 

To analyze the concurrence of the ground states we note that 
Hxx{h) and the z projection of the collective spin operator, 
S z = YJ)=i °/i are simultaneously diagonalizable, i.e. 



^JaW+i+Joyj+i+hS^St 



= 0. 



(23) 



This is not true for the anisotropic interaction. The fact that S z 
is a conserved quantity implies that while the eigenvalues of 
Hxx(h) are functions of «, the eigenstates are independent of 
h. We label the (possibly degenerate) ground state of Hxx{h) 
with its corresponding eigenvalues ^ of S z , as |\|4 ; (n)). Notice 
the expectation value of S z in terms of the fermionic number 
operators is 



(s z ) = (iLi(i-2«I^)) 

= n -2(r k=l Al k A 2k ) 
= n-2(r k= A<A2k) 



(24) 



Introducing the concurrence bilinear form C n (|\|/),|(|))) 
PW, satisfying C n {\\?)) 

Cn{S z \\ft{h)),S z \titW)) = 



= |C„(|\|/),|\|/))|, we find 

slc n {\Vi{h)),W s i{h))) 

{^i{h)\Sll5S z U- l UWi(h)) 

(Yi(h)\si(-s z )i3\Yi(h)) 



= -s 2 z c n (\yf(h))M(h))). 

(25) 
Here we have used the reality of the eigenstates and the eigen- 
values s z , and the fact that S z is time-reversal antisymmetric. 
Consequently, 



Cn(\y$(h)) 



ifs^O. 



(26) 



The n-concurrence of the ground state of Hxx (0) is one as 
per Prop. III. II therefore it is identified with the quantum num- 
ber s z = 0. As the magnetic field is increased, the energy gap 
decreases between the state |v|/5(ft)) and a state with differ- 
ent symmetry, |V/ C0)> unt il me Y become degenerate at a 
magnetic field strength n cr „. At this value of n, the ground 
state n-concurrence is zero. For \h\ > h C rti, it is energetically 
favorable for the spins to align meaning the spin projection 



\s z \ can only increase so the n-concurrence remains zero by 
Eq. |5S] In order to identify the degeneracy point we must 
consider two cases: the situtation with n/2 even, and n/2 odd. 
For both we assume J > 0, although a completely analogous 
argument can be made for J < 0. 



A. Case n/2 even 

At zero magnetic field and for J > the ground state energy 
is the lowest eigenvalue of H + : 



Eq = —8/ > cos — * 

k=i V " 

corresponding to the ground state 

n/2 

\^(\h\<h crit )}= n <-i^ ( 2H)i°)- 

A-=n/4+l 



-4/csc(7t/n) (27) 



(28) 



As h increases from zero, Eq becomes degenerate with the 
lowest eigenvalue of H : 

"/ 4 /n2k\ 
E min = -8/ ^ cos - 2/i -4J= -4/cot(7t/«) - 2h 

k=i \ n J 

(29) 
corresponding to the state 



W^ihcrit)) 



r "/2-1 

1 1 ^2*^-2* 

k=n/A- 1 



Al|0). (30) 



The magnetic field strength where this degeneracy occurs is 
h crit = — 2/[cot(7t/n) — csc(7t/n)] = 27 tan ( — j. (31) 



B. Case n/2 odd 

At zero magnetic field and for J > the ground state energy 
is the lowest eigenvalue of H : 

(»- 2 )/ 4 / % 2k\ 
Eq = -8/ £ cos = -4Jcsc(n/n) (32) 

fc=i V n J 

corresponding to the ground state 

re/2-1 



Wl{\h\ < h crit )) 



1 1 A 2lA-2k 

k=(n+2)/4 



Aj|0). (33) 



As h increases from zero, Eq becomes degenerate with the 
lowest eigenvalue of H + : 



(»-2)/4 

-87 y\ cos 

k=l 



%{2k-\) 



- 2/i = -4/cot(jc/n) - 2/i 
(34) 



corresponding to the state 

l¥ A 7 2 (^wr))= h A| t _ 1 Al (2Jt _ 1) |0>. (35) 

k=n/A+\ 

The magnetic field strength where this degeneracy occurs is 
again 

' 71 



hcrit = 2/ tan 



In 



(36) 



IV. THERMAL ENTANGLEMENT 

The thermal state of a system with intrinsic Hamiltonian H 
is given by 



pr = -j- = LL — i—^ 

^ j k=\ * 



(37) 



where p = 1/IcbT, with kg the Boltzmann constant, and the 
{|X^)} w& m e energy eigenstates of H corresponding to en- 
ergy Ej with degeneracy g } -,. Here Z = Tr[e~P H ] is the partition 
function. 

In the case that the intrinsic Hamiltonian is time reversal 
symmetric, the spin-flipped thermal state satisfies 

l5p T l5- 1 = I5e-^ H 13- 1 /Z 

= e-^ aau ~ 1 /Z (38) 

= Pr- 

Thus {A,;} = specKpTlipT^ 1 ) 1 ' 2 ] — spec[p;r] is just the set 
of probabilities to be in a thermal eigenstate. Using Eq|5]we 
then find 

C n (p T )=imxio,2e-^/Z-l\, (39) 

where Eq is the ground state energy of H. Because the par- 
tition function is a sum of positive terms, it is confirmed that 
C(pr) = if the ground state of H is degenerate. 

The Hamiltonian Hxx(h) is not time reversal symmetric 
for h^O, However, the concurrence for thermal states with 
a finite magnetic field is simply proportional to the concur- 
rence at zero field as we now show. Using the notation 
Z[h) = Tr[<rP* x W], we find 



{p T Up T U ) i - ( z(h) V z(h) U 



r^xxW e -po% x (*)«- 



1/2 



Z(h) 



Z(h) 



1/2 



,-WxxW e -Wxx(->>) 



1/2 



(40) 



-2P%x(0) 
Ah) 1 



Zjfi)PT\h=0, 



1/2 




C n { PT ) 




-S -6 -4 

ln(fe fl T/J) 



FIG. 2: Plot of the ^-concurrence of thermal states in 100 qubits 
of the Hamitonian Hxx{h) as function of magnetic field strength 
h. Shown are plots at several temperatures in units of the interac- 
tion strength: t = kgT/J. As the temperature drops to zero, the n- 
concurrence as a function of h approaches a step function with a dis- 
continuity at the critical magnetic field strength h cr i t . For n = 100, 
hcru = 27 tan(jt/ 100) « 0.03 147. 



where pr|/,=o denotes the thermal ensemble at h = 0. In de- 
riving Eq. |40|we have used Eq. [7]in the third line and the fact 
that [Hxx(0),hS z ] — in the fourth line. The concurrence is 

C„( P r) = |||max|o,2 e -P £ VZ(0)-l|. (41) 

where Eq is the ground state energy of Hxx (0). 

For the isotropic XX model, the partition function is calcu- 
lated using the the exact diagonalization in Sec. |HI]| 13] 



Z(h) 



2"- 



;nI=iCosh 2 (e 2 jt- 



IP)- 



n^=isinh 2 (8 2 ^ 



-n£=i ' cosh 2 (e 2<r P)cosh(e p)cosh(e„p) 



iP) 



„/2- 
11*= 1 



sinh 2 (e 2i p ) sinh (e p ) sinh (e„ P ) ) . 



(42) 

Using this expression for Z(/i) and the value of the ground 
state energy in Eq. [^ the concurrence is readily computed, 
(see Figs. |2]and[3]) 

It is useful to know how mixed the quantum state can be 
before quantum correlations are lost. For this purpose we de- 
fine a threshold temperature T t i, as the temperature where the 
concurrence changes from a positive quantity from below T,/, 
to zero above 7J/,. The value of p r /, = (ksT t h)^ 1 is obtained by 
solving the equation 



2e 



-P// ; £o 



/Z(0) -1=0. 



(43) 



An approximate solution for T t i, can be found in the limit of 
large n. The ground state energy Eq. 15 1 Iper particle converges 
to, 



lim — 



lim 



-47csc(7t/n) 



47 

K 



— (44) 



In the same limit, £ 2/ t_i « e 2 ^, and the logarithm of the parti- 
tion function per particle converges to 



lim 



lnZ(0) _ 2 f 71 / 2 



7t 7o 



atoln[2cosh(27Pcos(co))]. (45) 



FIG. 3: Plot of the ^-concurrence as a function of logarithm of the 
temperature for thermal states of the quantum XX model with zero 
magnetic field. The concurrence is plotted for several values of the 
total number of spins n. 



We can reexpress the integral as an infinite series in (7p) 



-l. 



lim 



lnZ(0) 



= 47p/7t- 

= 47P/JI- 

= 47P/JI- 

w 47P/JI- 

= 47P/7T. 



l^ l2 d&\n[\+e- AJ ^°<^] 



2(-if 



— (-i) t+1 



f* /2 d(de- 4kJ V C0 < a 



=1 ^7T 



Io(4fc7p)-L (4fc7p)] 

v ~ (-l) m + 1 r(l/2+m) 

1 nk ^™=0r(l/2-m)(2W|3) 2 »'+l 
V oo (-l) m (l-2- 2m - 1 )^(2m+2)r(l/2+m) 



jtT(l/2-m)(2/p) 



2/n+l 



(46) 

where Irj (x) (Lo (x) ) is the zeroth order Bessel(Struve) function 
and ^(x) is the Riemann Zeta function. The approximation in 
the fourth line is valid when 7p 3> 1 1 1711 . Inserting the two 
asymptotic expressions into Eq. |^] and neglecting terms that 
fall off like (7p)~ 3 and faster, we find 



hTti, 
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Ml 



(47) 



One can then conclude that entanglement as measured by the 
n-concurrence is greater than zero when the temperature is 
less than the interaction energy per particle. The derived ex- 
pression for T,h underestimates the exact value T,/, (exact) by 
roughly 20% (see Table|H) 

We compare our result with a related threshold temperature 
studied in 1 14] for a ring of qubits coupled via the antiferro- 
magnetic Heisenberg interaction. There it was found that the 
threshold temperature for nearest neighbor 2-concurrence ap- 
proaches a constant in the thermodynamic limit. This seems to 
indicate that quantum correlations between interacting qubit 
pairs are insensitive to changes in the character of the global 
thermal state as the system size increases. In contrast, for the 
XX model, the threshold temperature for the n-concurrence 
decreases inversely with the number of qubits. While we 
study a different spin chain Hamiltonian, qualitatively we ex- 
pect a different threshold temperature for n-concurrence ver- 
sus 2-concurrence in any spin chain model. This is because 
any time reversal symmetric Hamiltonian, of which the XYZ 
model is one type, has a maximally concurrent ground state 



TABLE I: Ratio of the approximate analytic value of threshold tem- 
perature T,/, from Eq. 1571 to the numerically computed exact value 
T p ,'° (exact) for periodic (p) and open (o) boundaries. The ratio is 
computed for several decades of the system size n. 



n 


10 


10^ 


10^ 


10 4 


%„ 


0.8159 


0.7960 


0.7959 


0.7959 


T/,: (exact) 


T,„ 


0.9979 


0.8969 


0.8887 


0.8879 


t,i (exact) 



for any even n provided it is non-degenerate. The size depen- 
dence of T t h is a consequence of the fact that the ground state 
population decreases with the size of the system for a fixed 
temperature. 



V. OPEN BOUNDARIES 

The previous analysis was done for the quantum XY model 
in one dimension with periodic boundaries. In most experi- 
mental situations it will be easier to construct the spin chains 
with open boundaries. In this section we rederive the critical 
magnetic field and threshold temperature for this case. 

The quantum XX model in one dimension with open 
boundaries is 



Hxx(h) = IfclJt&p} 



i-i. 



+ i 



■°5°/h)+*I3. 



i°} 



= Z%i2J(ajaJ +1 



■ojoj +l )+hZU^r 



(48) 

It is assumed that n is even. After the transformation from 
Pauli operators to fermionic operators, the Hamiltonian as- 
sumes the form, 



n-\ 



H X x (A) = nh - 27 £ (a j H 
;'=i 



■a.jaj+i) — 2hy\a-aj (49) 



7=1 



The Hamiltonian is a quadratic form on the vector space 
of the creation and annihilation operators {aj,a:}. Written 
as an n x n matrix in the basis {a\,a2,---a„}, Hxx(h) is a 
sum of two terms. The first is proportional to the identity 
and the second is a tridiagonal matrix with elements — 2h on 
the diagonal and elements —27 on the nearest off diagonal. 
The eigenvalues and eigenvectors of the tridiagonal can be 
solved for exactly yielding Hxx(h) — Y*k^kb\bk + C, where 
A^ = — 1h — 4Jcos(%k/(n + 1)), bj,bj are the quasi-particle 
creation and annihilation fermionic operators, and C is a con- 
stant. The value of the constant is determined by demanding 
that Tr[Hxx(h)] = 0. Summing over all possible particle oc- 
cupations, each particle energy At appears £"=n ("7 ) = 2" _1 
times and the constant satisfies — 2"~ L£=i A^ + 2"C = or 
C = nh. The eigenenergies can be expressed compactly as 



At zero magnetic field and for J > 0, the ground state energy 
is 



Eq = -4/e£cos(ji*/(b + 1)) 



= -47 cos 



W+T) j csc W+T) I sin 



(n+2)K 
4(«+l) 



(51) 
The corresponding ground state is an eigenstate of S z — 
Yfj=i ^i with eigenvalue s z = 0. As the magnetic field is in- 
creased from zero, this eigenstate becomes degenerate with 
the s z = —2 eigenstate at the critical value 



he 



27sin(7i/(2(« + l))). 



(52) 



For \h\ > hcrit, the n-concurrence of the ground state is equal 
to zero. Notice for n 3> 1, h cr u ~ Jn/n which is the same 
asymptotic as in the XX model with periodic boundaries. 

For the n-concurrence of thermal states, we need an ex- 
pression for the partition function Z(h) = e - ™"'). In the 
representation of the quasi-particles, the Hamiltonian can be 
written as a tensor product: 



^xx(«)=0 



k=\ 



The partition function is then 



z(h) = nu^"( 



$h(„-$A k . 



-P(A*H 




1) 







(53) 



e -P(^+LL 1 A t/2)ri^ =1 2cosh(pA jt /2) 
= 2 n nLi cosh (P(« + 27cos(rafc/(n + l))) 



(54) 



Given an expression for the partition function, the concur- 
rence can be calculated for any magnetic field and temperature 
using Eq. ^\ The threshold temperature can be computed by 
the same proceedure used in Sec. II VI In the large n limit, the 
ground state energy at zero magnetic field is 



lim — 



47 



(55) 



which is the same asymptotic as the case with periodic bound- 
aries. Again in the large n limit we find 



lim 



lnZ(0) _ 2 f n / 2 



n— >°o n % Jo 

and the threshold temperature is therefore 
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rfoorn[2cosh(2/Rcos(co))]. (56) 



ksTth 



nil 



(57) 



In Table|l]the exact value of the threshold temperature is com- 
puted numerically for comparison with the open boundaries 
case. 



VL MEASUREMENT 



£'- { Y* A k\ Afc€{-2ft-4/cos(jcfc/(n + l)),0} \+nh 

■ k=l J 

(50) 



For any measure of entanglement in a many body system 
it is important to find a practical method to construct its cor- 
responding observable or set of observables. Computing the 



n-concurrence on pure states is generically hard because the 
time reversal operator is not physical and therefore does not 
correspond to any single observable on the system. It has been 
shown that for pure states the n-tangle, which is the square of 
the n-concurrence, can be computed in terms of the multipho- 
ton Stokes scalar 11511 . There are an exponentially large num- 
ber of multiphoton Stokes parameters that need to be mea- 
sured in order to compute this scalar. Nevertheless, it may 
still be more efficient than resorting to direct state tomography 
over all n qubits which generically requires 4" measurements. 

A notable feature of Hamiltonians over spin chains such as 
the XX model is the reality of the Hamiltonian. This implies 
that the energy eigenstates are real and the concurrence of the 
ground state is then 



C„(|Y,» = |<Y*|U|Y«>| = \{V g \Xl(-i<?j)\Vg 

7=1 



(58) 



Measuring the expectation value of the many body oper- 
ator S = nn=i( — *'°/) can be done in various ways and the 
preferred technique will depend on the particular system. 
One approach would be to introduce an ancillary qubit pre- 
pared in the state \+ x )a = 1/v2(|0),a + | 1)a) and use this an- 
cilla as a control in a sequence of n controlled rotation gates 
rR = l Aaj(— iO y ). After measuring the ancilla in the \± x )a ba- 
sis, the concurrence is obtained from the measurement proba- 
bilities as |Prob(A == 1) -Prob(A ==0)| = C„(|\|/ g )). 

Measuring the n-concurrence for mixed states is much more 
difficult and will probably require some amount of state to- 
mography. In the case of thermal samples, if the temperature 
can be measured by some means then the results in this paper 
place a bound on how high it can be before quantum correla- 
tions as measured by C„ are lost. The analysis above shows 
that nearly maximal n-concurrence is achievable in the XX 
model with zero magnetic field for sufficiently small temper- 
atures. Because the temperature must be less than the interac- 
tion energy per particle this will be experimentally challeng- 
ing to observe. 

One candidate system may be cold trapped neutral atoms. 
Recently, it was proposed to simulate the quantum XY model 
in a lattice of spins trapped in an optical lattice llal . The lat- 
tice can be designed to trap an antiferromagnetic array of spin 
polarized atoms and the interactions can be engineered via 
"always on" ground state collisions between nearest neighbor 
atoms in the lattice. The advantage of using trapped atoms 
is that the lattice can be loaded from an atomic Bose Ein- 
stein Condensate (BEC) which are routinely prepared at ex- 
tremely low temperatures (~ 50 nK). For example, in a ID 
lattice with 100 qubits at a temperature of 50 nK, one would 
need interaction strengths of J m h2n x 16KHz in order to 
be in the regime of non-zero n-concurrence and an interac- 
tion strength J = ft2n x 115KHz to be close to maximal n- 
concurrence. Generating a sufficiently strong interaction be- 
tween neutral atoms is challenging, however, for appropriately 
tuned trapped potentials the collisional interactions can be on 
the order of a few tens of kilohertz. 



VII. CONCLUSIONS 



We have investigated the behavior of many qubit entangle- 
ment in one dimensional spin chains. The entanglement mea- 
sure used, the n-concurrence, is a global measure in the sense 
that its value goes to zero if any qubit is completely disentan- 
gled. This measure is appropriate in the context of spin chains 
because it probes the onset of a break in the time reversal sym- 
metry when a magnetic field is introduced. For pure states 
in n even qubits, the n-concurrence is maximal at zero mag- 
netic field. In the XX model, the entanglement is not immedi- 
ately lost when a magnetic field perturbation is added, rather 
it jumps abruptly to zero at a critical magnetic field strength 
corresponding to the first degeneracy point. Using the sym- 
metry properties of the Hamiltonian under conjugation by the 
time reversal operation we are able to obtain expressions for 
the entanglement of thermal systems of arbitrary size n. The 
entanglement was shown to be finite below a certain critical 
magnetic field strength and threshold temperature. 

Several outstanding issues remain. First, the objection may 
be raised that maximal n-concurrence in the ground states of 
Hxy (0) does not imply global entanglement whatsoever be- 
cause subglobal entanglement, i.e. tensor products of time 
reversal symmetric states, also have maximal C„ (see Sec. |IIJ. 
However, it has been proven JUlI mat the ground state of the 
XX model has non vanishing 2-concurrence over all pair sep- 
arations. Hence, there can be no way to partition the ground 
state into disentangled subsets of qubits and therefore the en- 
tanglement is truly global. Second, it is of interest to study 
how inhomogeneous couplings affect the n-concurrence. Per- 
haps the behavior of global entanglement could give a signa- 
ture of local defects. Finally, because time reversal is not a 
spatial symmetry, it may be possible to investigate entangle- 
ment in the usually more complicated model of spin chains 
embedded in d > 1 dimensions. 

Another question for future research is whether the highly 
entangled ground state of a spin chain can be efficiently trans- 
formed into a state useful from the point of view of quantum 
information processing (QIP). Spin chains are apromising ar- 
chitecture to implement many QIP tasks such as entanglement 
distribution 1 18], quantum state swapping 1 19], and quantum 
computation |20]. Preparing and observing many -body en- 
tanglement will be an important first step to realizing these 
challenging tasks. It is appealing to try to generate useful 
entanglement by allowing the system to naturally cool to its 
ground state. The results obtained in this paper help define 
limits to the environmental conditions under which one type 
of entanglement (the n-concurrence) can be prepared in this 
way. It is known |9] that all states of fixed n-concurrence are 
orbits of the Lie (symmetry) group K = SO(2") of all unitary 
evolutions which admit time reversal antisymmetric Hamil- 
tonains. In particular, there is some (and there are many) 
Hamiltonian(s) H with H = -UHU~ l and |GHZ) = k\\?g) 
for k = exp( —iHt/h ). It is unclear whether there exist effi- 
cient operations in K to transform from one to the other. 
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APPENDIX A: CONCURRENCE OF MIXED STATES 

This appendix expands upon published work of Wootters 
and Uhlmann Il7l l23ll . recalling their closed-form expression 
for the minimum Eq. [5] It is, to our knowledge, the first 
treatment of this problem that is complete and self-contained. 



1. Notation and conventions 

For either scalars or vectors, we denote complex conjugates 
by an overline. Also, we often use both real and complex 
operators. Thus, we forego bra-ket notation in favor of the set 
{ej} jIq of standard basis vectors, i.e. column matrices with 

(ej)t = Sj. The symbols W ,x i and C pxc > denote pxq real, 
complex matrices respectively. Finally, we use a dagger for 
adjoint throughout. For example, e 1 : is formally equivalent to 
a bra if the operator is intended to be complex. 



Then it is standard that Sp(£,R) n 0{2£) = U(£) |EI pg. 447, 
Lemma2.1(c) of Chap.X.2] as we have verified above. 

We also note generalities on time-reversal symmetry op- 
erators, denoted 0, on a Hilbert state-space of complex di- 
mension N. Such a is an M.-linear map of the state 
space possessing (i) complex anti-linearity, (ii) orthogonality 
r = I2N, and (iii) projectively involutivity © 2 = (e'^I^R, 
<j) ^ 0. Given that x is the R-linear map corresponding to 
complex conjugation on C , we note that any antilinear map 
(hence 0) may be written as = C0rX for co € C NxN . Due 
to orthogonality of and x, COr is orthogonal. Hence CO is 
unitary, given C-linearity. We further claim that e I(p is ±1. 
Indeed, 2 = co r xcorX = corCOr, so that coco = e l,f I N . But 
Tr(coco) = Tr(coco), demanding that the trace is real. Hence 
e l<f is ±1. We refer the +1 case as bosonic time-reversal sym- 
metry operators, and the — 1 case are fermionic time-reversal 
symmetry operators. 



2. Background 



3. Analysis of 



£&■"*; 



Throughout, we distinguish C and R as vector spaces 
over R. Thus, throughout we use following choice of R- 
isomorphism: 



If w 



x = a + ib € ' 



1,21 



(Al) 



iv € C is a C-linear map, we note the corre- 



sponding R-linear map under the isomorphism. Specifically, 

w (a + ib) = (u + iv) (a + ib) = (ua — vb) + i(va + ub) for 



-iv £ 



>txl 



S R 2ix2e (A2) 



Finally, we note that xx is not xt^x^. Rather, by Pr we intend 
R-linear map that extends the C-linear map given by p. Com- 



puting, as x 
ib T ) = (aa T ■ 



■ ib £ C , we see p = xx 1 = (a + ib)(a 



-bb') + i( 



ab ), for 



The following proposition will be quite important to what 
follows. Thus, we present a careful argument. It improves an 
earlier inductive proof of the second author and was suggested 
by Dianne O'Leary. 

Proposition A. 1 Let Xq > \\ > X2 > ■ ■ ■ > Ki-i > be an 
ordered set of real numbers. 



ifLj = \ X j <x 0' 



then 



e-i r i-\ 

A-o — V Xj = min < ) , e' 'X 



;6cvfy_ie: 



(A5) 



IfYj ; = 1 Xj > A,o, then there exist {Qj} . =0 so that 
i-\ 



(A6) 



Y T V 



-ive 




(A3) 



We also should breifly comment on the structure of mr 
for uu^ = U, u €E C (x( . Clearly any unitary map of C^ will 
lift to an orthogonal map of R , but this is not the com- 
plete structure. Complex multiplication by i lifts as a matrix 

T „ ] . Then an orthogonal map o of M. 2e has o = mr 

for some unitary u iff oJ = Jo. Equivalently, a unitary map is 
an orthogonal map which is also C-linear. Adding a tad more 
language, as one speaks of the unitary group U(£) there is also 
a real symplectic group [21, pg. 446] given by 



/ = 



Sp(£, 



{Ae 



p 2^x2^ . 



A T JA 



J} 



(A4) 



The proof uses the followig two lemmas. The first is a stan- 
dard result in combinatorics. 

Lemma A.2 Let Xq > X\ > X2 > ■ ■ ■ > %i-\ > 0, and label 
S = {0, 1 , • • • , £ — 1}. Then there exists a partition S = S1US2 
such that 



*o > E ^ - E ^ > 

jeSi J£S 2 



(A7) 



Proof: Let S = Si U S2 be a partioning chosen to minimize 



max< T^jeSi ^j — EjeS 2 X j^(- Label 5 > to be this minimal 

difference, i.e. 5 = ZjeSj h - I.jes 2 ^j- 

Assume by way of contradiction that 8 > Xq. In particular, 
some element of 5i must be nonzero, say Xk > 0. Now place 
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§i=Si- {X k } and S 2 = S 2 U {A*}. Then label 8 > by 






\8-2X k 



(A8) 



Now either |8 — 2Xk\ — 8 — 2X k or else |8 — 2Xk\ = 2X k — 8. In 
the former case, note that 8 — 8 = 2A# > 0, i.e. 8 > 8 contra- 
diction. In the latter case, 8 — 8 = 28 — 2X k > 2(8 — Xq) > 0, 
contradiction. □ 

Lemma A.3 Let to>t\> 0, and /afoe/ L(0) = |fo + e l9 *i|. 
Then for every s G [fo — fi,fo + ^i]» f^ere is *ome 0o € [0,7c] 
smc/z f/zaf L(8o) = s. 

Proof: Clearly L(0) is continuous. Now L(7t) = fo — f i an d 
L(0) = to + ti, so that the result follows by the Intermediate 
Value Theorem. □ 

Proof of Prop. IA.ll We begin with the first item. Any 



{0j r =1 must have 



Ejii&Xj 



<lfjZ\\e m iXj\.Thxis 



l-i 



> le^Aol- 



E^A ; - > Xo-Z X J (A ^ 



On the other hand, L,=i Xj < Xq demands the last expression 
of Equation lA9l is contained within the set being minimized. 

For the second part, label the (£ — lj-index set S = 
{1, • • • ,£ — 1}. By Lemma |A~21 we may partition S = S1US2 
so that 



A > Ai > £ Xj <- £ Xj > 



Label fo = HjeSi 'Xj an ^ ? i 



£/es 2 V 



(A10) 

Now fo + fi = 



£*•_] A-j > Aq. Hence by Lemma lAJl there is some so that 



|?o + e' e fi| 



I>, + £e% 
jeSi jes 2 



= Ao 



Hence for some e'^, we have e'^Xo = YijeS-i Xj + TijeS 2 e ' 9 A;'- 
Thus e ! >Ao + I /G5l e'%- + £jes 2 e ;(e+,l) A, =0. D 



Remark A.5 As always, we could produce a left action 
by considering the right action of the adjoint. However, we 
find that confusing in this context. Checking the right action, 
consider {wy}j~J = (£ • m 1 ) • u 2 for m 1 , u 2 unitary. Then 



t-it-i 



l-\ / N-\ 



W J = E Y, X P U )k U kl = H x p[ H u )k u ll ( A12 ) 

k=0p=0 p=0 V k=0 



This is £ ■ (u l u 2 ). We also remark that if £ is an ensemble for 
p, then so likewise is £ ■ u. Indeed, 



L%lwjw) = I}loltoIp=o^^^4 Mt ^4 



(-1W-I 



I ;= oI i= 0^ e ; MM e 



' jyj/ T -Trl — 



i A k 



L,k=0 X k x k ~ P 



(A13) 



Thus, the U(£) action respects the density matrix structure. 
It has been proven [22] the action is transitive on the set 
returning p; any subnormalized ensemble for p of length £ 
arises in this way. Every p possesses an ensemble of length 
n, due to the spectral theorem. We also remark that since 
dimcEndc(C Ar ) = N 2 , it is in some sense wasteful to take 
£ > N 2 . However, the arguments would not simplify with this 
convention. Finally, note that evidently any ensemble for p 
must have length at least rank(p). 



(All) Definition A.6 Let £ = {xj}j = q be any ensemble of a fixed 



density matrix p. The preconcurrence c 2p 



c 2p( 



x'&x = X' COX, 1^0 



0, x = 



(A14) 



4. Concurrence of Ensembles 

We set further notations regarding density matrices. These 
make later arguments more clear and terse. 

Definition A.4 Let p e C NxN be a density matrix, i.e. p = 
P > P > 0, and Tr(p) = 1. A subnormalized ensemble for p 
of length £ is any collection of vectors £ = {*/}tJ C C^ 
satisfying p = LjZo x j x f This is associated to a normalized 
ensemble of length £ given by {ymZj, where yj — if Xj = 
and yt — Xj/\xj\ else. We further define a right action of the 
unitary group U(£) on the set of subnormalized ensembles of 
length £ as follows. If u — (uj k ) eC ext , uu< — 1(, then £ • u = 

{wjYjZl where wj = YlZ$ x ke)ue k = Y, e k Z l x k u jk . 



Notethatforx^O, |c 2p (x)| = |x| 2 C 2/ ,(x/|x|). We then define 
the concurrence of an ensemble £ to be the following sum, 
for £ = {**}££: 

G,(£) = £!«,(%)! = E \A 2 Ci P {x/\x\) (A15) 

We also make the following definition, fixing some preferred 
ensemble £0 for p. 

C|_(p) = min{C2,;(£ ) ; £ is an ensemble of length I for p} 



= min{ C 2p (£o ■«);«£ U(£) } 
We then define C 2p (p) = min{ CUp) ; £ > 1 }. 



(A16) 
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Theorem A.7 (| 231) View a given time-reversal symmetry op- 
erator G ^-NxiN ^ Endm.{pin)> and denote the M.-linear 



map corresponding to p as Pr £ 



x2N 



. Since is antiuni- 



tary, in particular orthogonal, & T = . We write = COrT 
for CO G i^NxN an £ ^ tne complex-conjugation viewed within 
Endg,(C N ). Moreover say {hj} I J=Q = spec[M(p)] are the 
(real) eigenvalues of 

M (p) = (\/P«pco + VP) 1/2 > notingthat 
[(VptopcotVp) 1 /2] R = (V p i 0p K 0r V p i) i/2 ^') 



ordered so that A-o > A-i > X 2 • 
and bosonic, we must have 



>X N -i >0. Then if N = 2" 



C 2p (p) = maxl OM - £ Xj | (A18) 



In view of Proposition lA.il we prove this theorem in two 
steps. Put T — max{0,A-o — Y?j=\ **]}• ^ n tne f" rst ste P> we 
produce an ensemble £ mm for p (of length AO such that the 
T = C2p(£ mm ). In the second step, we restrict to case r^O 
and show that any ensemble £ for p has C 2p (r, )>T. Each of 
the two steps is organized within a subsection, culiminating in 
Pronositions lA. 14llAT5l and lA.171 



5. Existence of a Minimizing Ensemble 



Lemma A.8 Let A £ 

oo = I(, so that (oAo 



sl'xf 



MM : 



Then there is some o G 



vtxt 



{oAo 1 



(oAo )t-\i-\. 



1 yit-1 



Sketch: Let T = \ Y,j=o a jj> tne average value of the diagonal 
elements of A. Since the trace of oAo T coincides with the 
original trace, we seek to set all diagonal elements to be T . If 
every diagonal element is T already, then take o = Ig, 

Else some diagonal element exceeds T, and some is less 
than T '. By choosing an appropriate permutation matrix n 
and relabeling 1X411^ as A, we may suppose without loss 
of generality (WLOG) that am < T, a u > T. Now la- 
bel R(t) — e~ ! '° ®Ii—2< Consider the continuous function 
t h^ [R(t)AR(t) T ]oo = cos 2 (f)aoo - sin 2 (f)an. The Intermedi- 
ate Value Theorem shows that for some fo the diagonal entry 
is T . Now induct on the number of entries equal to T . □ 

Lemma A.9 Let r| £ C (x(: be a complex, perhaps non- 
Hermitian matrix: r| = T| . Then there exists a unitary matrix 
u G U(£) so that 

ut\u t = A, A = diag(ko,Xi,...,Xg-i) (A19) 

Moreover, we may choose Xq > Xi > • • • > A-£_ i > 0. 

Remark A.10 This is not a diagonalization of r|! Indeed, 
u is unitary; u~ l = u' ^ u T . The point of this proof is to 
diagonalize the M-linear map t|rT, which is symmetric. <) 
Proof: Since T| is symmetric, we must have the following: 



TlR = 



T T 

U = U , V = V 



(A20) 



Now if fj{ denotes scalar multplication by i, then clearly /j,r| = 
t|u ; . Thus 7riR = r| K 7, for J = (jh)m = {-io y )®k. 

Now let fj be the associated complex anti-linear map, i.e. 
f[x = V[x. Then t|r = T|rX, and noting that X = diag(/Ar, —In) 
produces 



(f|)if 



riRi = 



(A21) 



Hence fjn = fj^, i.e. f|R is diagonalized by some matrix 
orthogonal matrix. We next argue that such an orthogonal 
matrix o may be chosen to be symplectic (within Sp(£,M).) 
Per earlier discussion, o £ Sp(£,M) n 0(2£) will then demand 
o = Mr for some u G U(£). 

Step # 1: Note that for any eigenvalue of fJR, say X, we also 
have —X as an eigenvalue. For Jt\m = — tJr7, since fj is C- 
antilinear. Hence T[gJv = — /tJrv = (— X)Jv, given v £ V\ C 
R 2e . 

Step #2: Choose a collection of positive eigenvalues and 
take A = diag(A-o, . . . ,Xt-i). For the corresponding orthonor- 
mal eigenvectors vo, ■ ■ • , V£- 1 , consider the orthogonal matrix 
o = (vo- • • v^_i7vo' • 'Jvi-l)- Then7o = oJ, i.e. o £ Sp(£,R). 
Moreover of\ino T = A (—A). Note that without loss of gen- 
erality, A reflects a choice Xq>X\ >■■■ X(-\ > 0. 
Step #3: Let u be the unitary matrix resulting from o in the 
last step. Then the final equation demands u^u* = At, i.e. 
uyyiu^x = A, i.e. ur\u T = A. □ 

Definition A.ll Let % = {xj\ -Z^ be an ensemble for some 
density matrix p, and let © = COT be a time-reversal symmetry 
operator. Then we define r| (■£ , 0) to be that matrix whose en- 
tries are [t|(e , &)]jk = xMxk = XyCOx^. We often suppress the 
arguments and write r| = T| (E , 0) when the context is clear. 

Lemma A.12 We have the following basic properties of 
T1(E,0). 

1. r|(£ •«,©) = mt|(e,©) u*. 

2. For any two ensembles "E\, TL*i of a given p of length £, 
we have 

spec[r\('Ei,®)r\(£i,&)] = specify,®) f\(£2,&) ] 

(A22) 

3- If£ = {xjYjll, then c 2p (x k ) = r\ kk . 
Proof: Keeping £ = {xj}*f~Q, note that £ • u is Wj = 

Hk=0 X l< e ] ue k- Thus 



[T|(E -u,®)] jk = w]aw k = 
lfp=llfgllel^ejxl(ax g elue g = 

i±oEio%ft(£.©)W"i 



(A23) 



P<l u qk 



The first item results, comparing entry by entry. For the sec- 
ond item, recall that £2 = El • u for some unitary matrix u. 
The third item is the definition of preconcurrence. □ 

Proposition A.13 Fix a density matrix p G C f x£ . Suppose 
is bosonic, so that CO = C0 r . 
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1. Then there exists a subnormalized ensemble Eo = 
{xj} -Zq so thatr\{'Eo,&) is diagonal, real, and contains 
only nonnegative entries. 

2. Label M(p) = [^/pcopco 1 ^] 1 / 2 . If the length of a given 
ensemble £ is N, then 

spec[M{p)} = spec[y/T\(-E,®)T\(-E,®)] (A24) 

Proof: We first prove ItemlO Fix £ of length £ and 0, so that 
T| <G Cr x . Since CO = C0 r , we have r| = r| r . Hence, there exists 
by Lemma IX. 91 a unitary matrix u such that u T|(£ ,©) u T — 
A, A = diag(A,o,A-i,--- ,Xt-\), and Ao > Ai > ••• > Xi-\ > 0. 
Thus,A = r|(£ -m,0). 

For Item|21 let £ = {*/} wj ■ Now by Lemma lA.121 we may 
suppose without loss of generality that £ is a subnormalized 
eigenensemble, so that jtjjtjt = Sy and pxj = VjXj for {Vy}^T 
the set of eigenvalues of p. Form the unitary matrix w with 
wej = Xj, i.e. w — L,=o x j e ]- Then 



wr)riw T = ny^xje) JZ^epy/Vjxlafm 

= t L^o y/VpXpxl) ^ ( ^=o >ffi*A 
= v/pcopco+Vp 

(A25) 
In the last line, we note that {xj}j_ is a normalized ensemble 

for p. Hence, spec( rjff ) = spec( A /pcopco t v / p). D 

Proposition A.14 Let Xq > Xi > ••• > A,jv— i > be an order- 
ing of the spectrum ofM(p) per Theorem \A.7\ Suppose in ad- 
dition that T = A-o — £/=i Xj > 0. Then there exists a subnor- 
malized ensemble £ OT / n of length N such that C2 P { r E m i n ) — T. 



Proof: By Proposition IA. 131 we produce an ensemble £o 
of length N so that r|(£o,@) = A = diag(Ao,Ai,- •• ,Xn-i). 
Moreover, {A,/}^7q = spec[M(p)], since rjff = A 2 . Now 
label the phase matrix <t> = diag(l, ;',;', • • • ,i), and note 
that T = Tr[Or|(£o,©)<I> + ] = Tr[r|(£o ■ *,©)]■ Moreover, 
<Px\('Eq, ©)<!>' is real. Hence, by Lemma lA~8l there exists an 
orthogonal matrix o £ C NxN such that [or|(£o ■ &,@)o T ]jj = 
Tl(Eo • (®o),®)jj = T/N, < j < N - 1. We claim that 
we may now take £ mm = £o ■ ( < J > °)- Indeed, put £ mm = 
{yj}p0- Then y]® yj = r\('E min ,®) JJ = T/N for each < 
j <N- 1, so that 



C2 P (£ mi n) = Ep, b;fC2 P (j;/b;|) 
E-L"o bj^i =Z N ^T/N= T 

This concludes the proof. 



(A26) 

D 



Proposition A.15 Lef AT = 2", anc/ /abe/ Ao > Ai > • • • > 
Aat_i > an ordering of the spectrum of M(p) per Theorem 
IA. 71 Suppose in addition that T = Xq — Y,;=i ^j < 0- Then 
there exists an ensemble T, nl ^ n for p such that Cip{'E m i n ) = 0. 



Proof: Again, Pronosition lA. 1 3l produces an ensemble £o of 
lengths so thatr)(£o,®) = A = diag(Ao,Ai, ■ ■ ■ ,^/v-i) with 
{XjYI_}— spec[M(p)]. In this event, we appeal to Proposi- 
tion [AT] to assert that there must exist a collection of angles 

~JXj\. WLOG, take 9 = 0. 



iJV-l 



-JV-l„rt 



Now put 4> 

Trfl^Eo,®)"^] = 0- Recall the Hadamard computation 

H = AsYri k?=o(~iy k \k){j\> which is unitary. Now consider 



such that = | Zj^q & 

diag(l,e-'' e '/ 2 , 



^-M N -i/2y We then have 



Tr[//»"4ni(Eo,©)^ T ^"] 



We next claim that £_ 



Tr[Ti(£ -(3># 8 "),e)] = 
(A27) 
mm = Eo • (*//® n ) = {?)}*£ is an en- 
semble for p consisting of concurrence zero states. Indeed, 

put {Zj} N jZo — £ o ■ *. Then z]®Zk = e^'Xfi). Now consider 
that due to the application of the Hadamard computation 



yj = E z kZj 

k=0 

for each £& = ±1. Hence, C2n(v/) = 0. 



(A28) 

□ 



6. Minimality of max{0, A l!J~i Xj } 

Lemma A.16 Let r = rank(p) > rank[M(p)]. Consider the 
first r, concievably nonzero eigenvalues Xq > Ai > ■ ■ • > X r > 
of M{p). Then for any ensemble E of p of arbitray length 
I, for T| = T|(e , ©), the r largest eigenvalues of ^/r|rf are also 
A >Ai >••• >X r -\ >0. 

Sketch: To begin, write E as {xj}jZq. Then Tiff = 
^^ x]copco t L^ Xi : Now note that I > r. Then by transi- 
tivity of the U(£) action, we may suppose without loss of gen- 
erality that E is a subnormalized eigenensemble, perhaps with 
trailing zero eigenvectors. Letting {v;}t be the normalized 
eigenensemble, again consider the matrix w wej — yj. Then 
again for {Vj}y =0 the nonzero eigenvalues of p, 



wT|r|w' 



-e-i. 



FjZoyje})r-J e P V^yl<<>p~rf 

N /ptopco t A /p 



(A29) 



Although w is no longer square, this statement nonetheless 
implies the result on truncated spectra. □ 

Proposition A.17 Let E be any ensemble for a density ma- 
trix p of length L Let T — max{ 0, Ao — £/=i Xj} for the Xj 
nonincreasing and coinciding with the spectrum of M(p) — 
(VpcopV^/p) 1 / 2 . Then C 2p {'E) > T. 

Proof: Fix r) = r|(£,@). We recall by Lemma lA9l that 
there exists some unitary matrix u so that ur\u T = A. By 
Lemma lA.161 we have A = diag(Ao,Ai,A2,- • • ,X(-\) for Ao > 
Ai > ■ ■ ■ > Xt-i > the £ largest eigenvalues of M(p) — 
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(y/pcopco^y/p) 1 ' 2 . Moreover every nonzero eigenvalue of 
M(p) appears within the set of the t largest. 

We first consider the case that T — Xq ~ E/=i^j > 
0. Now note that by the Definition IA.11I for T|, we 
have C 2p (£) = EJlJh^l ^'— ■• "■ -<..a..^ 

The result then 



Y!„2§Uj P 'kpUkp, so that r| 



77 



Moreover, T)^ = (uAu 7 )jk 
- L p =o u jp^p- 



follows from the Schwarz inequality, given Ll=o \ u 



Clp^E) = Iffj 

_ yC-l 



177 



r £-l v £'-l 2 J 
L p= oLj=i u jp^p 



> V 



Ep=0^-P 

L p =Ap 



Y e ~ l u 2 



1: 



(A30) 



This concludes the proof for A-o > L,=i Xj. 



Thus, suppose T = 0. Then the statement is vacuous, since 
always C2p("E ) > 0. This concludes the proof. □ 
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